A Multi-Code Analysis Toolkit for Astrophysical Simulation Data

The analysis of complex multiphysics astrophysical simulations presents a unique and rapidly growing set of challenges: reproducibility, parallelization, and vast increases in data size and complexity chief among them. In order to meet these challenges, and in order to open up new avenues for collaboration between users of multiple simulation platforms, we present yt (available at http://yt.enzotools.org/), an open source, community-developed astrophysical analysis and visualization toolkit. Analysis and visualization with yt are oriented around physically relevant quantities rather than quantities native to astrophysical simulation codes. While originally designed for handling Enzo's structure adaptive mesh refinement (AMR) data, yt has been extended to work with several different simulation methods and simulation codes including Orion, RAMSES, and FLASH. We report on its methods for reading, handling, and visualizing data, including projections, multivariate volume rendering, multi-dimensional histograms, halo finding, light cone generation and topologically-connected isocontour identification. Furthermore, we discuss the underlying algorithms yt uses for processing and visualizing data, and its mechanisms for parallelization of analysis tasks.Comment: 18 pages, 6 figures, emulateapj format. Resubmitted to Astrophysical Journal Supplement Series with revisions from referee. yt can be found at http://yt.enzotools.org

Effects of Varying the Three-Body Molecular Hydrogen Formation Rate in Primordial Star Formation

The transformation of atomic hydrogen to molecular hydrogen through three-body reactions is a crucial stage in the collapse of primordial, metal-free halos, where the first generation of stars (Population III stars) in the Universe are formed. However, in the published literature, the rate coefficient for this reaction is uncertain by nearly an order of magnitude. We report on the results of both adaptive mesh refinement (AMR) and smoothed particle hydrodynamics (SPH) simulations of the collapse of metal-free halos as a function of the value of this rate coefficient. For each simulation method, we have simulated a single halo three times, using three different values of the rate coefficient. We find that while variation between halo realizations may be greater than that caused by the three-body rate coefficient being used, both the accretion physics onto Population III protostars as well as the long-term stability of the disk and any potential fragmentation may depend strongly on this rate coefficient.Comment: 29 pages, 7 figures; Accepted for publication in The Astrophysical Journa

AGN Clustering in the Local Universe: An Unbiased Picture from Swift-BAT

We present the clustering measurement of hard X-ray selected AGN in the local Universe. We used a sample of 199 sources spectroscopically confirmed detected by Swift-BAT in its 15-55 keV all-sky survey. We measured the real space projected auto-correlation function and detected a signal significant on projected scales lower than 200 Mpc/h. We measured a correlation length of r{sub 0} = 5.56{sup +0.49}{sub -0.43} Mpc/h and a slope {gamma} = 1.64{sup -0.08}{sub -0.07}. We also measured the auto-correlation function of Tyep I and Type II AGN and found higher correlation length for Type I AGN. We have a marginal evidence of luminosity dependent clustering of AGN, as we detected a larger correlation length of luminous AGN than that of low luminosity sources. The corresponding typical host DM halo masses of Swift-BAT are {approx} log(M{sub DMH) {approx} 12-14 h{sup -1}M/M{sub {circle_dot}} which is the typical mass of a galaxy group. We estimated that the local AGN population has a typical lifetime {tau}{sub AGN} {approx}0.7 Gyr, it is powered by SMBH with mass M{sub BH} {approx}1-10x10{sup 8} M{sub {circle_dot}} and accreting with very low efficiency, log({epsilon}){approx}-2.0>. We also conclude that local AGN galaxies are typically red-massive galaxies with stellar mass of the order 2-80x10{sup 10} h{sup -1}M{sub {circle_dot}}. We compared our results with clustering predictions of merger-driven AGN triggering models and found a good agreement